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ABSTRACT 


This  study  explores  the  dynamics  of  diffusive 
convection  which  occurs  when  cold,  fresh  water  overlies  warm 
and  salty  water.  Diffusive  convection  is  generally  observed 
in  high-latitude  regions,  particularly  in  the  Arctic  Ocean, 
between  the  base  of  the  mixed  layer  and  the  top  of  the 
Atlantic  water.  The  primary  convective  regime  in  this 
region  is  characterized  by  the  spontaneous  formation  of  well 
mixed  layers  separated  by  thin  high-gradient  interfaces 
known  as  thermohaline  staircases. 

Data  analysis  and  analytical  considerations  are  used  to 
estimate  the  vertical  heat/salt  mixing  rates  and  their 
dependencies  on  the  large-scale  environmental  parameters. 
The  ice-tethered  profiler  data  from  the  Beaufort  Gyre  are 
analyzed  to  determine  the  origin  of  the  thermohaline 
staircases  and  the  mechanism  for  selection  of  the  preferred 
layer  thickness.  Based  on  the  analysis  of  Beaufort 
staircases,  we  suggest  that  the  layer  thickness,  as  well  as 
the  vertical  heat/salt  fluxes,  is  controlled  by  the  patterns 
of  merging  events  in  which  relatively  small  steps  are 
systematically  eliminated.  Significant  concerns  are  raised 
with  regard  to  the  direct  extrapolation  of  laboratory 
derived  flux  laws  to  ocean  conditions  —  a  conventional 
approach  to  the  analysis  of  the  Arctic  staircases.  An 
alternative  and  more  accurate  method  of  analysis  is  proposed 
which  involves  recalibration  of  the  laboratory-derived  flux 
laws  for  the  oceanic  conditions. 

Extrapolated  diffusive  convective  fluxes  are  in  the 
range  of  1-6  WmA-2,  which  is  comparable  to  magnitude  of 


v 


fluxes  currently  unaccounted  for  in  the  Arctic  heat  budget. 
We  propose  that  the  parameterizations  of  the  diffusive 
fluxes  in  thermohaline  staircases,  developed  in  this  thesis, 
can  be  used  to  enhance  understanding  of  Arctic  climate 
changes  and  predictive  capabilities  of  large-scale  numerical 
models.  Preliminary  findings  are  indicative  of  the 
importance  of  diffusive  convection  for  sound  propagation  in 
the  Arctic  region  —  the  problem  of  great  interest  for 
various  Naval  research  applications  in  the  area. 
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I. 


INTRODUCTION 


The  polar  ice  cap,  consistently  regarded  as  a  hindrance 
to  data  collection  in  the  Arctic,  is  currently  under¬ 
sampled.  With  evidence  of  increasing  Arctic  temperatures 
and  decreasing  sea-ice  thickness,  the  need  to  parameterize 
unresolved  ocean  processes,  and  include  them  in  climate 
models,  is  becoming  a  necessity.  This  study  focuses  on  the 
double  diffusive  region  of  the  Beaufort  Gyre. 

A.  DOUBLE  DIFFUSION 

Before  discussing  the  Arctic  data,  the  basic  concepts 
behind  double  diffusion  should  be  reviewed.  Although  double 
diffusion  occurs  in  many  different  systems,  it  has  a  vast 
impact  on  ocean  dynamics.  Turner  (1973)  provided  an  entire 
chapter  on  the  dynamics  of  double  diffusive  regimes  in  his 
book  on  the  buoyancy  effects  in  fluids.  He  describes 
double-diffusion  as  the  process  related  to  the  presence  of 
more  than  one  component  contributing  to  the  density  of  the 
ocean  with  different  molecular  dif fusivities .  In  the  ocean 
case,  the  dominant  density  components  are  temperature  and 
salinity.  Ocean  double  diffusion  takes  two  forms.  The  most 
studied  is  the  salt  finger  regime,  where  warm,  salty  water 
lies  above  cold,  fresh  water.  The  other  regime,  which  is 
the  focus  of  this  study,  is  diffusive  convection  where  cold, 
fresh  water  overlies  warm,  salty  water.  Stratification 
patterns  which  support  salt-fingering  and  diffusive 
convection  are  shown  in  the  schematic  diagrams  in  Figure  1 . 
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Figure  1.  Simplified  diagram  of  ocean  double  diffusive 

regimes,  a)  Salt  Finger  Regime  b)  Diffusive  Convection 

Regime . 


1 .  Salt  Fingers 

The  majority  of  research  has  focused  on  the  salt¬ 
fingering  regime.  Stommel  et  al .  (1956)  brought  interest  to 
the  salt  finger  regime  by  hypothesizing  that  by  exploiting 
the  warm,  salty  water  over  cooler,  fresher  water,  an  energy 
source  termed  "the  perpetual  salt  fountain"  could  be  formed. 
The  concept  being  that  a  thin  membrane  in  the  shape  of  a 
tube  could  allow  temperature  diffusion,  but  not  salinity 
diffusion,  and  create  an  energy  source  by  making  a  parcel  of 
salty  water  rise  as  it  tries  to  reach  a  state  of 
equilibrium.  While  the  proposed  idea  has  yet  to  be  fully 
developed  in  the  ocean,  the  interest  in  salt  fingering  made 
its  debut. 

Stern  (1960)  described  salt-fingers  as  instabilities 
that  arise  when  relatively  warm,  salty  water  lies  over 
colder,  fresher  water  of  a  higher  density  and  consists  of 
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long  fingers  of  rising  and  sinking  water.  The  warmer,  salty 
water  rapidly  loses  heat  while  retaining  its  salinity.  The 
increasing  parcel  density  makes  it  sink  further,  drawing 
more  warm,  salty  water  from  above,  resulting  in  the  fingers 
of  salty  fluid.  Figure  2  shows  the  resulting  salt  fingers. 
In  this  case  the  magenta  dye  is  applied  to  the  warm,  salty 
water . 


Figure  2.  Laboratory  simulation  of  salt  fingers.  The 

magenta  plumes  are  salty  water  pushing  into  the  colder 
fresher  water.  (From  University  of  Alaska  Fairbanks 
website,  www.ims.uaf.edu). 


2 .  Diffusive  Convection 

The  dynamics  of  double-diffusive  convection  differ 
substantially  from  those  of  salt  fingering.  Diffusive 


convection  occurs  in  both  a  smooth  gradient  region  and  a 
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stepped  configuration,  which  is  more  common  in  nature 
(Figure  3)  .  For  the  diffusive  convection  case,  relatively 
cold,  fresh  water  overlies  warm,  salty  water  with  a  thin 
interface  in  between. 


Cold,  Fresh 

Thin  interface 


Ah 


Warm,  Salty 


Figure  3 . 


Basic  set-up  of  diffusive  layers. 


Molecular  fluxes  across  the  interface  are: 


ccT 'tm  K.j 


a  AT 
A  h 


and  pFSM  k. 


pAS 

Ah 


where  a  is  the  coefficient 


of  volume  expansion  due  to  unit  temperature  change,  P  is 
the  proportional  density  change  produced  by  unit  density 
salinity  change,  FTM  and  FSM  are  the  molecular  fluxes  of  heat 
and  salinity  respectively,  and  kt  and  ks  are  the  molecular 
dif fusivities  of  heat  and  salt,  respectively. 


Since  kt^>ks,  for  aAT  and  PAS  of  the  same  order  of 
magnitude,  the  density  flux  Fp  =  pPFSM  -  paFTM  <  0.  The  rapid 

molecule  diffusion  of  temperature,  relative  to  salt, 
produces  a  downward  density  flux  across  the  interface.  The 
region  immediately  below  (above)  the  interface  becomes 
denser  (lighter) ,  which  maintains  the  top-heavy  convection 
in  mixed  layers.  Figure  4  illustrates  the  physical 
mechanism  for  diffusive  layering. 
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Figure  4.  Physical  mechanism  behind  diffusive  layering. 


There  are  several  regions  of  the  world' s  oceans  with 
the  potential  for  double-diffusive  convection,  including  the 
Arctic  Ocean  (Kelley  et  al .  ,  2003)  .  The  intensity  of 
diffusive  convection  is  controlled  by  the  so-called  density 
ratio,  which,  for  the  diffusive  convection  case,  is  defined 
as  : 


,  = 

Lp  a  AT 


(1) 


where : 


1  (dp 
p{dT 

/?=.! lJd^ 


/  p* 


p 


\dS' 

J  p*  p* 


(2) 


Fluxes  decrease  with  increasing  density  ratio.  The 
range  l<i?^<10  is  most  susceptible  to  diffusive  convection 
(Kelley  et  al .  ,  2003)  .  Analysis  of  Arctic  density  ratios 
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shows  that  the  majority  of  the  Arctic  Ocean,  specifically 
the  Beaufort  Gyre,  has  density  ratios  in  the  3-10  range. 
(Figure  5) 


Figure  5.  Regions  susceptible  to  double  diffusion.  The 

light  grey  areas  indicate  where  density  ratio  is  from 
3-10.  Dark  grey  indicates  l<i?/9<3.  The  red  circle 
indicates  the  Beaufort  Gyre. 


The  Arctic  is  diffusively  stratified  due  to  the  process 
of  ice  formation,  brine  rejection,  and  subsequent  melting. 
The  process  creates  the  cold,  fresh  layer  that  lies  on  top 
of  the  warm,  salty  water  carried  from  the  Atlantic  Ocean  to 
the  Arctic,  through  typical  ocean  circulation.  Figure  6 
illustrates  Arctic  currents  and  circulation  patterns. 
Figure  7  provides  a  vertical  cross  section  of  the  Arctic 
Ocean . 
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Figure 

(i 


6.  Arctic  circulation  pattern,  resulting  in 

conditions  favorable  for  diffusive  convection, 
llustration  by  Jack  Cook,  Woods  Hole  Oceanographic 

Institution . ) 
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Figure  7.  Cross  section  of  Arctic  water  mass.  Cold, 

fresh  water  overlies  warm,  salty  Atlantic  water. 

(Illustration  by  Jayne  Doucette,  Woods  Hole 
Oceanographic  Institute . ) 

In  Figure  8  Fernando  (1987)  simulated  the  layering 
process  by  heating  a  salinity  gradient  from  below.  The 
diffusive  layers  are  indicated  in  the  diagram.  Similar 
laboratory  simulations  have  also  been  conducted  using  a 
sugar  solution  as  the  denser  liquid.  Figure  9  is  a  more 
completely  evolved  set  of  diffusive  layers  in  the 
laboratory.  Turner  (1968)  used  this  shadowgraph  to  show  the 
scale  of  layers. 
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Figure  8.  Side  view  of  diffusive  layering  created  by 

heating  a  salt  gradient  from  below,  after  Figure  6 
(Kelley  2003)  A  is  the  first  quasi-stationary 
interface.  B  is  the  secondary  interface,  forming  above 

it . 


Figure  9.  A  series  of  convecting  layers  and  heat-salt 

diffusive  interfaces  (Turner  1968)  . 
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B. 


THEORETICAL  BACKGROUND 


Turner  (1965)  hypothesized  that  if  steps  are 
sufficiently  large  (H— »oo),  the  dependence  of  the  vertical 
fluxes  of  temperature  and  salinity  (  FT  and  Fs)  on  H  would 
be  weak,  and  therefore  fluxes  would  be  determined  by  the 
temperature  and  salinity  variations  (ATLAS')  across  the  step: 


Fj  Fj(aAT ,  pAS  ,Kj  ,v,  g)  FT(aAT  ,Rp,KT,g,z,Pr) 


Fs  =  Fs  (aAT,  PAS,  kt,ks,  u,  g)  =  FT  (< aAT ,  Rp,tcT,  g,  r,  P ) 


(3) 


where  u  is  viscosity,  r  is  the  diffusivity  ratio  and  Pr  is 
the  Prandtl  number 


'o' 


\KT  J 


Using  the  dimensional  argument  that  a  non-dimensional 
number  can  depend  only  on  other  non-dimensional  numbers. 
Turner  then  suggested  that  the  so-called  Nusselt  number. 


Nu  = 


Ft 

AT 


measuring  the  fluxes  would  be 


related  to  the 


Raleigh  number,  measuring  boundary  forcing: 


gaATTP 

kto 


Nu  =  C(Ra)n  =>  aFT  = 


r  gaATH3  Y 
V  ktv  J 


'  kt 


aAT 

H 


From  (3) ,  we  know  that  there  should  be 


(4) 

no  dependence  on 


H  in  (4),  therefore  n=l/3  or  aFT  ~  (aAT )3  .  These  equations 
are  called  the  "4/3  Flux  Laws".  (Turner  1973) 
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Another  non-dimensional  quantity  of  interest  is  the 


flux  ratio: 


Y  = 


aFT 


which  measures 


potential  energy  lost  by  temperature 
goes  into  the  potential  energy  of  salt 


the  fraction  of 

stratification 
stratification . 


the 

that 


Turner  (1965)  ,  argued,  on  the  basis  of  dimensional 
considerations  and  laboratory  experimentation,  that  for 
given  molecular  properties,  the  flux  ratio  is  a  function  of 
Rp  .  Experiments  yield  a  striking  result  that  for  Rp>2  the 
flux  ratio  remains  constant  to  within  the  experimental 


error . 


In  the  laboratory  this  mean  value  is 


iL. 

aFT 


0.15  ,  with 


maximum  heat  density  flux  for  a  given  AT  occurring  at  Rp=2  . 

Figure  10  is  a  plot  buoyancy  flux  ratio  versus  density  ratio 
in  the  laboratory. 
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Figure  10.  Buoyancy  flux  ratio  RF(Rp)  for  diffusive 

convection  in  the  laboratory.  Data  from  three 
different  laboratory  sources  Crapper  (1975),  Turner 
(1965),  and  Newell  (1984)  from  Kelley  (1990).  The 
solid  line  is  the  empirical  fit. 

The  aforementioned  models  offer  simple  conceptual 
explanations  for  the  dependences  of  vertical  fluxes  on 

temperature  and  salinity  variations  across  the  steps 

(ATLAS')  .  It  follows  that  knowledge  of  what  determines  AT 
in  the  staircases  is  necessary.  Kelley  (1988)  and  Radko 

(2005)  suggest  that  a  layer-merging  process  controls  step 
characteristics.  This  merging  process,  which  will  be 
discussed  further  in  Chapter  III,  can  be  defined  as  the 

coalescence  of  thin,  adjacent  layers  to  form  new,  thicker 
layers,  and  their  eventual  equilibration.  The  major  goal  of 
this  study  is  to  examine  the  validity  of  flux  law  models  and 
predictions  of  merging  theory.  For  that,  we  turn  to  Arctic 
data . 
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C.  ICE-TETHERED  PROFILER 

Originally,  thermohaline  staircases  were  not  visible  in 
ocean  profiles.  This  was  due  to  data  taken  only  at  standard 
levels  versus  a  continuous  profile  of  the  region.  In  order 
to  characterize  the  diffusive  convective  layers,  high- 
resolution  data  are  necessary.  This  study  uses  high 
resolution  data  collected  from  the  Woods  Hole  Ice-Tethered 
Profiler  Project. 

In  conjunction  with  the  Beaufort  Gyre  Exploration 
Project,  a  team  at  Woods  Hole  Oceanographic  Institute  (WHOI) 
has  launched  several  Ice-Tethered  Profilers  (ITP's)  in  the 
Beaufort  Sea.  Analysis  of  diffusive  convection  primarily 
uses  data  from  ITP1  and  ITP  3. 

ITP's  are  moored  profilers  that  are  designed  to  make 
daily  measurements  of  water,  pressure,  salinity 
(conductivity) ,  and  temperature  in  a  region  where  ice 
impedes  access  and  collection  of  this  data.  Designed  to 
last  for  three  years  when  deployed,  as  long  as  the  ice  floe 
survives,  the  profiler  sends  its  measurements  and  GPS 
location  to  a  surface  float,  which  is  then  relayed  back  to 
WHOI.  The  profiler  is  made  up  of  a  small  surface  capsule,  an 
800  meter  line,  a  moored  profiler,  and  a  small  anchor.  The 
profiler  is  modified  so  to  fit  through  a  10"  hole  in  the 
surface  (Krishfield  et  al . ,  2006). 
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Figure  11.  Schematic  of  Ice-Tethered  Profiler 

(Krishfield  et  al .  ,  2006) . 

The  high  resolution  (25  centimeter)  data  posts  to  the 
WHOI  ITP  website  for  further  analysis  in  its  raw  form. 
There  are  future  plans  to  de-spike,  filter,  and  correct  the 
data  using  a  sensor  response  correction  (Johnson,  2007)  . 
These  corrections  should  occur  in  the  fall  of  2007  and  be 
considered  for  future  work. 

ITP  1  was  the  second  of  seven  ITP's  deployed.  It  was 
launched  in  the  Beaufort  Sea  on  August  15,  2005  at  78°  51.1 
N,  150°  15.9  W  during  the  JWACS  2005  cruise  on  the  CCGS 
Louis  S.  St.  Laurent.  ITP1  was  deployed  on  a  4.6  meter  thick 
multiyear  ice  floe  in  a  location  designed  to  drift  through 
the  Beaufort  Gyre  Observing  System  (BGOS)  array  over  the 
following  year.  ITP3,  on  August  23,  2005  at  77°  36.1  N, 
142°  11.8  W,  during  the  same  cruise.  It  was  deployed  on  a 
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3.5  meter  thick  multiyear  ice  floe  in  a  location  designed  to 
drift  through  the  Beaufort  Gyre  Freshwater  Experiment  (BGFE) 
array  during  the  following  year.  Both  units  operate  on  a 
fast  sample  schedule  of  4  one-way  profiles  between  10  and 
760  meters  depth  each  day. 


ITP  drift  (yellow  line)  and  latest  location  (triangle), 

BGOS  moorings  (white  circles)  and  annual  ice  drift  from 
IABP  (grey  vectors)  on  IBCAO  bathymetry  (shading). 

Figure  12.  ITP  1  and  ITP  3  drift  tracks,  (a)  Drift  track 

for  ITP  1  and  (b)  is  ITP3  drift  track.  From  Woods  Hole 
Oceanographic  Institute  Ice  Tethered  Profiler  site 
(http: //www.whoi . edu/itp/data . html ) . 
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II.  CHARACTERISTICS  OF  THERMOHALINE  STAIRCASES 


A.  STEADY  CHARACTERISTICS 

The  temperature,  pressure,  and  conductivity  raw  data 
were  used  to  create  temperature-depth  and  salinity-depth 
profiles  through  the  diffusive  layer,  characterized  by 
prominent  staircases.  In  the  Beaufort  Gyre  Region, 
"staircase"  structures  were  prominent  from  depths  of  about 
200  to  400  meters.  Therefore,  all  analyses  were  computed 
within  that  depth  range.  Figures  13  and  14  are  temperature- 
depth  and  salinity-depth  profiles  illustrating  the  staircase 
characteristics.  As  denoted  in  Figure  13,  H  is  defined  as 
the  distance  between  the  midpoints  of  adjacent  interfaces. 
Figure  14  indicates  AS  measurements.  AT  is  found  in  the 
same  manner.  All  vertical  calculations  were  conducted  with 
this  data.  The  step  is  defined  as  the  vertical  rise  between 
two  adjacent  interfaces.  a  and  jB  were  calculated  using 
the  MATLAB  seawater  toolbox.  Although  they  vary  with  depth, 
average  values  of  a  are  6*10A-5  K_1  and  jB  are  8*10A-4 

(PSU)-1.  Average  AT  was  0.04  degrees  Celsius,  and  average 

(  jBAS^ 

AS  values  were  0.01  PSU.  The  density  ratios  = -  , 

{  p  aATJ 

computed  for  the  ITP  data,  ranged  from  3-7.  Table  1  shows 
the  data  for  ITP  1  and  ITP  3  averaged  by  month. 
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Pitfle1153 


Figure  13. 


Temperature  Depth  Profile.  The  section  of 
the  profile  is  highlighted.  H  calculations  are 
determined  from  the  midpoint  of  each  staircase  depth. 
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Profiler  1153 


Salinity  (psu) 


Figure  14.  Salinity  depth  profile.  Same  location  as 

previous  figure.  (Salinity  is  calculated  using  MATLAB 
polyfit  to  de-spike  data) . 
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AT 

H(m) 

Rp 

ITP1 

September 

2005 

0.0403 

4.1635 

3.2713 

October 

0.0421 

4.7896 

3.6998 

November 

0.0437 

4.6067 

3.8795 

December 

0.0466 

4.2825 

4.759 

January 

2006 

0.0437 

4.1593 

3.8973 

February 

0.0434 

3.8276 

3.8536 

March 

0.0456 

4.3723 

3.7838 

April 

0.0406 

4.0872 

3.8512 

May 

0.0397 

3.4183 

4.0875 

June 

0.0429 

3.9029 

4.1846 

July 

0.0472 

4.31 

4.3198 

August 

0.0441 

3.9876 

4.6595 

September 

0.0435 

4.0714 

4.9589 

October 

0.0466 

4.2825 

4.759 

ITP3 

September 

2005 

0.0373 

3.2653 

5.5181 

October 

0.0432 

3.6059 

4.548 

November 

0.0394 

3.5775 

4.6802 

December 

0.0394 

3.4346 

4.8329 

January 

2006 

0.0378 

3.0995 

4.7914 

February 

0.04115 

3.3285 

4.638 

March 

0.0455 

3.8186 

4.4491 

April 

0.046 

3.5688 

4.2637 

May 

0.0455 

3.6922 

4.3631 

June 

0.0477 

4.1085 

4.734 

July 

0.0356 

2.7624 

5.4397 

August 

0.0391 

3.1684 

5.3817 

Average  over  both  profilers 

0.0426 

3.8342 

4.446 

Table  1 . 


Average  data 
over 


from  ITP  1 
the  entire 


and  ITP  3  by  month  and 
data  set . 
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Figure  15  plots  aAT  (non-dimensional  temperature 
variation  with  depth)  as  a  function  of  H  (calculated  step 
height.)  Although  only  the  month  of  January  is  presented, 
similar  results  were  obtained  for  each  month.  From  this 
comparison,  it  is  theorized  that  the  change  in  temperature 
is  relatively  linear,  and  is  nearly  independent  of  the 
height  of  the  step.  The  slightly  negative  slope  is  caused 
by  the  dependence  by  a  on  pressure.  In  general,  H  values 
increase  from  about  2  meters  at  a  depth  of  240  meters  to 
above  20  meters  at  a  depth  of  400  meters,  while  AT  remains 
on  average  of  0.04  degrees  C.  This  observation  is 
consistent  with  Turner's  suggestion  that  the  vertical  fluxes 
of  heat  and  salt  are  independent  of  step  height  (Equation 
(3)  )  .  Indeed,  if  the  vertical  fluxes  FT  and  Fs  are 
vertically  homogeneous,  then  AT,  AS  =  constant.  Figure  16 
shows  a  similar  relationship  between  salinity  and  H. 

x  io'6  January  2006  (Potential  Temp  ITP1) 


Figure  15.  Non  dimensional  aAT  is  plotted  as  a  function 

of  layer  thickness,  H  (meters) . 
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Figure  16. 


Non-dimensional  J3AS  is 
of  layer  thickness  H 


plotted  as 
(meters) . 
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In  Figures  17  and  18,  aAT  and  J3AS  are  presented  as  a 


function  of  density  ratio,  Rp  . 


Figure  17.  Comparison  of  aAT  and  Rp  using  average  data 

from  September  2005  to  October  2006. 
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x  10  s 


Figure  18.  Comparison  of  J3 AS  and  Rp  using  average  data 

from  September  2005  to  October  2006. 


In  Figure  19,  the  density  ratio  over  the  interface  is 
compared  with  the  density  ratio  over  the  step.  For  this 
comparison,  density  ratio  is  first  calculated  over  the  step, 
and  then  calculated  using  two  points  on  the  adjacent 
interface.  The  density  ratios  for  the  interface  and  steps 
are  then  averaged  by  month  and  plotted.  The  increased  value 
of  Rpinterface  over  RpStep  is  an  indication  of  the  ambiguity  in  the 

definition  of  density  ratio,  and  leads  to  a  possible  source 
of  error  in  the  4/3  flux  laws. 
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ITP1  data  averaged  by  month  from  September  2005  to  October  2006 


Figure  19.  Comparison  of  density  ratios  over  the 

interface  and  density  ratios  over  the  step. 


There  is  general  agreement  in  the  literature  that  the 
double-diffusive  temperature  and  salinity  fluxes  are 
controlled  by  the  density  ratio.  While  this  relationship  is 
known  to  exist,  there  is  still  considerable  uncertainty  in 
both  the  dynamics  and  values  of  fluxes  in  the  thermohaline 
staircases  using  a  straightforward  extrapolation  of  the  lab- 
based  flux  estimation.  Padman  and  Dillon  (1987)  calculated 
thermohaline  fluxes  in  the  Beaufort  Gyre  to  be  about  0 . 1 
WmA-2  for  2<Rp<l .  A  comparison  between  the  thermohaline 

staircases  that  Padman  and  Dillon  explored  and  the  current 
data  indicate  that  in  the  20  years  between  data  collection, 
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there  were  significant  changes  in  the  water  mass.  This 
leads  to  changes  in  the  calculated  fluxes. 

In  1987,  the  depth  range  of  the  staircases  was  300-430 
meters  with  A7"  on  order  of  0.01  degrees  C  vice  the  0.04 
found  today.  There  is  a  more  fundamental  difficulty  with 
regards  to  the  earlier  estimates  of  fluxes  in  diffusive 
staircases.  In  addition  to  the  differences  in  water  column 
characteristics,  Padman  and  Dillon's  calculations  were  made 
using  direct  extrapolation  from  laboratory  numbers.  In  the 
salt  finger  regime,  such  an  extrapolation  was  long  ago 
proven  unreliable.  In  the  diffusive  convection  regime, 
Kelley  (1990)  also  suggested  that  direct  application  of  the 
4/3  Flux  Law  to  ocean  conditions  is  equally  unreliable.  In 
order  to  validate  this  theory,  Radko  (2007) 's  merging 
theorem  was  applied  to  two  components,  temperature  and 
salinity.  These  theoretical  values,  using  laboratory 
results,  were  compared  to  the  actual  data.  The  result  was 
confirmation  that  the  4/3  Flux  laws  cannot  be  directly 
extrapolated  from  the  laboratory  setup  to  diffusive 
convection  in  the  ocean. 
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B. 


EVOLUTIONARY  PATTERN 


Initially,  the  possibility  of  mergers  existing  in  the 
ITP  data  came  from  a  plot  of  conductivity  versus  time. 
Conductivity  was  used  vice  salinity  because  it  could  be 
taken  in  its  raw  form.  In  Figure  20,  lines  of  constant 
conductivity,  approximating  the  observations,  are 
superimposed  on  the  measured  data.  The  merging  events 
occurred  twice  in  this  figure  —  that  is  where  the  lines  of 
constant  conductivity  meet  and  create  fewer  lines. 


np  i 


200  400  600  800  1000  1200  1400 

Time 


Figure  2  0 . 


Examples  of  merging  events,  as  seen  from  the 
inspection  of  the  conductivity  data. 
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Overall,  we  identified  ten  clear  mergers  to  use  in  the 
following  (Chapter  III)  theory-based  analysis.  Several  of 
these  mergers  are  presented  in  Figure  21.  Visual  inspection 
of  all  merging  events  suggest  that  layers  merge  when 
slightly  stronger,  in  terms  of  temperature  and  salinity 
variations,  interfaces  grow,  while  weak  interfaces  decay  and 
ultimately  disappear.  This  is  the  evolutionary  pattern 
classified  by  Radko  (2007)  as  the  B-merger  scenario. 


Figure  21 . 
point 


Temperature  variation  in  time.  Each  data 
represents  temperature  on  a  particular  profile. 
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III.  MERGING  THEOREM  FOR  A  TWO  COMPONENT  FLUID 


Radko  (2007)  theorized  that  the  formation  of  diffusive 
layers  starts  with  multiple,  small-scale,  homogenous  density 
layers  that  merge  until  they  reach  equilibration.  There  are 
two  types  of  merging  processes.  Radko  defines  them  as  H- 
merger  and  B-merger.  The  H-merger  is  defined  when  slightly 
thicker  layers  grow  while  thin  layers  shrink  and  ultimately 
disappear  (Figure  22) .  In  the  case  of  B-merger  (Figure  23), 
stronger  interfaces  grow  at  the  expense  of  weaker 
interfaces.  While  theoretically  both  types  are  possible, 
data  indicates  that  only  the  B-merger  exists  in  the  ocean, 
therefore  the  included  derivations  are  limited  to  the  B- 
merger . 


H-Merger 


t=3 


Figure  22 . 


Schematic  of  H-merger  evolution  in  time 
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B-Merger 


t=2 


t=3 


Figure  23. 


Schematic  of  B-merger  evolution  in  time 


The  derivation  for  one  component  and  one-dimensional 
mergers  in  Radko  (2007)  was  modified  to  two  components, 

salinity  and  temperature,  yielding  the  following  description 

of  the  B-merger  events.  This  is  accomplished  by  conducting 
a  linear  stability  analysis  for  a  series  of  identical  steps 
and  using  it  to  formulate  general  criteria  for  B-merging 
events . 

In  order  to  establish  the  basic  principals  of  the 
staircases  evolution  in  time,  we  consider  the  stepped 
solution  of  the  one-dimensional  conservation  equations: 

dT  _  dF 'T_ ' 

dt  dz  ,n, 

r  (5) 

dS  _  dF_  s_ 

dt  dz  _ 


The  stepped  solution  is  shown  in  Figure  24.  Figure  24a 
represents  a  staircase  consisting  of  identical  steps,  the 
basic  steady-state.  Figure  24b  represents  a  slightly 
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perturbed  state  in  which  steps  are  no  longer  equal.  The 
temperature-salinity  variations  across  the  two  steps  in 
Figure  24(b)  are  expressed  as: 


tu-t01=t-s 

T23  -Tn=T  +  S 
Su-Sm=S-s 


1S23  —  Sl2  —  S  +  £,  (<5,  £) « (T,  S ) 


(6) 


where  T  -  (dT /  dz)H  and  S  -  (dS  /  dz)H  are  the  variations  over  one 
step  in  the  undisturbed  staircase  in  Figure  24(a)  and  (t ,5) 
are  small  perturbations.  Assume  that  H  is  constant  and 
equal  where  H  is  the  distance  between  the  interfaces  at 
z  =  and  z  =  zH+1  . 
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Figure  24.  Schematic  diagram  illustrating  the  stability 

analysis  for  an  infinite  series  of  interfaces.  (a) 
Basic  state  consisting  of  identical  steps,  (b) 
Perturbed  state  in  which  the  T,S  jumps  at  the  even 
interfaces  are  slightly  decreased,  and  the  jumps  at  odd 
interfaces  are  increased. 


Then  integrating  (5)  over  the  interval  [zbot,  ztop]  where 
ztQp  and  Zbot  are  two  levels  which  do  not  vary  in  time: 


zbot  zbot 
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Using  the  identity: 


d 

ztop 

r 

7  qt  , 

[  Tdz 

=  f  — dz 

dt 

J 

dt 

zbot 

zbot 

d 

ztop 

r  _  , 

y  ds , 

Sdz 

=  f  — dz 

dt 

J 

_J  dt 

zbot  zbot 


Equation  (7)  is  rewritten  as 

J  ztop 

-jTdz=FT(z:J-FT(z,„) 

l  ( 9 ) 

r  ztop 

-j  Sdz=Fs(z,J-Fs(zt„) 

zbot  ^ 

Applying  the  integral  relation  in  equation  (9)  to 
individual  layers  [zn,zn+l]  yields,  where  Fj=F(Zi): 


J  zn+ 1 

zn 

7  Zn+ 1 

Zn 


The  major  contribution  to  the  integral  of  buoyancy  in 
(10)  comes  from  the  interior  part  of  layers  where  the 
buoyancy  gradient  is  nearly  uniform,  and  therefore, 


'«+ 1 

|  Tdz*T„tlH 

zn 
zn+ 1 

J  Sdz»S.MH 


(11)  where  Tnn+1  and  Snn+1  are  the 


temperature  and  salinity  values  at  the  layer  centers. 
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Applying  (8)  and  (9)  to  two  consecutive  layers  [zn,zn+ J 


where  n=l,  2  results  in: 

jt(nrn)=F; -f; 

tShs*)=f>-f~s 

jt(HSn)  =  F?-F? 

Using  equations  (6)  and  (12),  the  equations  for 
evolution  of  5  and  e  can  be  formulated,  and  then  simplified 
using  periodicity  conditions. 


H—  =  -2(F2T-FlT) 
at 

H  -j-  =  -2( F2s  -  Fts ) 
at 


(13) 


Considering  weak  perturbations  such  that  (S  T,s  <SC  H)  and 
linearizing  (15)  yields: 


(14) 
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At  this  point,  it  becomes  necessary  to  specify  how  the 
fluxes  in  (14)  depend  on  temperature-salinity  jumps  across 
the  steps.  For  that,  we  replace  the  instantaneous  fluxes  in 
interfaces  Fn  in  (15)  with  the  steady  'one-step'  flux 

F  (Radko  2007)  based  on  the  parameters  of  a  region  extending 
between  the  centers  of  two  adjacent  layers  (zn_ln  andzml)  : 


Fy T 

f2t 

Ff 

f2s 


*Ft(T12-T01,S12 

(^23  —  ^12  5^23 

*Fs(Tu-Tm,S,2 
(^12  _  ^01 » *^12 


^01’Z12 
*^12  ’  Z23 
^0 1  ’  Zl  2 
^01’  Z12 


-zoi) 

-Zl2) 

-zoi) 

_zoi)y 


(15) 


The  difference  between  fluxes  at  the  adjacent 
interfaces  in  (14)  reduces  to: 


8Ft 

.  8Ft  1 

Fj  -Ft  =25 

~  +  2s  ~ 

Z  1 

8T 

8S 

8FS 

.  8FS 

F?s  -F.s  =28 

~  +  2  £  ~ 

z  I 

8T 

8S  } 

and  (14)  becomes: 

d  _  4S_  dFr  4s  8F^_ ' 
~dt  ~  ~H~df  ~H~df 
d_  4  5  8F_ s_  4s  8F^_ 
dt£~  H  df  H  8S  , 


Substitution  of  the  normal  modes  (<?,£•,/?)  =  (<50,£-0,/z0)exp(/lf) 

in  (17)  yields  the  eigenvalue  equation  for  the  B-growth 
rates : 


+- 


H 


8Ft  8F 

+ 


8T  8S 


Ab  + 


16 
H 2 


8Ft  8Fs 
8T  8S 


8FS  8Ft  " 
df  8S  , 


(18) 
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If  the  free  coefficient  is  small  (to  be  expected  since 
it  is  proportional  to  the  flux  ratio  variation)  the  positive 
root  can  be  further  approximated  by: 

dP ’s  dFr  dFr  dF s_ 

x  4 w^ww  (19) 

H  dF^  dF^_ 

~W  +  ~dS~ 

The  tildes  denote  the  equilibrium  one-step  parameters. 
The  latter  affords  an  attractive  opportunity  to  validate  the 
merging  theory  (Huppert,  1971;  Kelley,  1988;  Radko,  2005)  by 
comparing  the  growth  rates  of  the  observed  merging  events 
with  the  theoretical  prediction  for  F(RP)  . 
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IV.  APPLICATIONS  OF  THE  MERGING  THEOREM  TO  OCEANIC 

DATA 


In  order  to  apply  the  merging  theorem  (19)  to  actual 
field  data,  explicit  expressions  for  fluxes  are  needed. 
While  little  is  known  about  the  flux  laws  in  the  ocean,  we 
assume  that  the  general  structure  of  the  laboratory-derived 
flux  laws  pertain  to  the  ocean,  but  additional  calibration 
of  amplitude  is  necessary.  Adopting  the  derived  flux  laws, 
but  recalibrating  them  for  oceanic  conditions,  (e.g.  Kelley 
1990)  yields: 


ccFr  =  AC 


gkj 


v  v  j 


J3FS  =  aFTy ;  y  = 


(aAT)4'3;  C  =  0.0032  exp(4.8/i?“ 72 ) 


Rp+\A{Rp-\f2 
\  +  \4(Rp-l)3'2 


(20) 


"A"  is  an  adjustable  coefficient  included  to 
recalibrate  Kelley-Turner  flux  laws  for  the  ocean. 


An  alternative  lab-based  formulation  was  suggested  by 
Marmorino  and  Caldwell  (1976) : 


HSP  =  (0.085)jfcr 


r  V/3 

ga_ 


AT 


4/3 


(21) 


H 


^—  =  0.101exp|4.6exp|^-0.54(i?p  -l) 

SP 


(22) 


Finally,  substitute  (20)  and  (22)  into  (19)  and  obtain 
specific  solutions.  Simplifying  (19)  takes  the  form: 


(23) 
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Now,  we  apply  the  ITP  data.  The  exponential  variation 
of  AT, AS  across  the  merging  layers  confirms  the 
interpretation  of  the  merging  as  an  instability  process. 
The  growth  rate  of  each  merger  is  determined  from  the  best 
fit  of  the  exponent  exp(/U)  to  (  ATa—ATaq  ),  where  ATa  ,  is  the 
temperature  variation  in  the  merging  layer,  and  A TA0  is  its 
value  at  the  beginning  of  the  merging  event.  Two  examples  of 
growth  rate  determination  are  provided.  In  the  first 
example,  the  density  ratio  is  slightly  higher,  leading  to  a 
smaller  growth  rate  and  a  longer  period  for  merging  (Figure 
25)  .  In  the  second  example,  the  density  ratio  is  only 
slightly  smaller,  but  the  growth  rate  is  twice  as  large, 
yielding  a  much  shorter  merging  period  (Figure  26) . 


123456789  10 

Time  in  Days 


Figure  25.  Change  in  temperature  vs.  time  during  the 

merging  period.  The  normalized  growth  rate  for  this 
merger  is  4.7E-04  and  the  density  ratio  is  2.8631. 
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Figure  26.  Change  in  temperature  vs.  during  the  merging 

period.  The  normalized  growth  rate  for  this  merger  is 
8.6E-04  and  the  density  ratio  is  2.4035. 


The  merging  growth  rates  are  then  normalized  as 
follows : 


A 


AH 


'normal 


(24) 


a  =  6e-5K~l 


(aAT)m 

where  H  is  the  initial  step  height  and  6(  =  De_;,A-  #  They  are 

then  compared  with  the  normalized  growth  rates  from  the 
theoretical/laboratory  values.  Kelley  (1990)  suggested  that 
the  coefficient  C  in  the  4/3  flux  law  differed  based  on 
density  ratio.  Based  on  laboratory  estimates. 


C  =  .0032  exp 


r  4.8  A 


nO.72 

\RP 


( 25 ) 


This  differs 
suggested 


from  Marmorino  and  Caldwell 

C=0. 085*0. 101*exp^4.6exp|^(-0.54|Rp-l|  j  . 
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(1976)  who 


(26) 


The  differences  in  the  calculations  of  laboratory  data 
lead  to  questions  regarding  the  precision  of  laboratory 
measurements.  In  order  to  estimate  the  error  of  the  model 
predictions,  both  Kelley  and  Marmorino  and  Caldwell's 
formulations  were  independently  used.  Figure  27  is  based  on 
Kelley's  expression,  which  required  a  coefficient  A=30.008 
to  fit  Beaufort  Gyre  data.  Marmorino  and  Caldwell's 
coefficient,  was  A=15.75  and  the  corresponding  prediction 
for  Anorm  is  shown  in  Figure  28.  Table  2  shows  a  numerical 
comparison  of  the  fit. 


Data  vs.  Theory 


Figure  27.  Using  the  coefficient  A=30,  theoretical  data 

fits  Beaufort  Gyre  data  within  95%  prediction  bounds. 
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x  ip'3  Fit  line  of  data  vs.  merging  theory 


Figure  28.  Using  the  coefficient  A=15.75  applied  to 

Marmorino  and  Caldwell  (1976),  theoretical  data  fits 
Beaufort  Gyre  data  within  95%  prediction  bounds. 


Rp 

AT 

H 

A  normalized 

A  normalized 
A=30 

A  normalized 
A=15.75 

1.9106 

0.1572 

24 

2.70E-03 

2.58E-03 

2.48E-03 

2.2509 

0.1007 

7 

3.1  IE-04 

3.1  IE-04 

7.93E-04 

2.2737 

0.087 

3 

3.60E-04 

3.60E-04 

7.43E-04 

2.4031 

9.61  E-02 

4.6 

1 .06E-03 

7.17E-04 

5.24E-04 

2.4035 

1.12E-01 

3 

8.63E-04 

7.17E-04 

5.23E-04 

2.4267 

1.18E-01 

5 

6.08E-04 

6.83E-04 

4.93E-04 

2.8631 

6.99E-02 

4 

4.71  E-04 

3.16E-04 

1.90E-04 

3.2997 

8.20E-02 

3.7 

1 .28E-04 

1 .75E-04 

9.25E-05 

3.3728 

7.91  E-02 

2.8 

2.81  E-04 

1 .60E-04 

8.35E-05 

3.3879 

9.09E-02 

3 

1 .90E-04 

1 .58E-04 

8.18E-05 

Table  2.  Data  comparison  of  calculated  normalized  growth 
rate  from  actual  data  and  the  normalized  growth  rates 
using  a  constant  factor  with  the  4/3  flux  law. 
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The  coefficient  A,  which  we  believe  is  necessary  to 
convert  from  laboratory  simulations  to  ocean  data  leads  to 

relatively  high  heat  fluxes  (HP)  for  low  density  ratios 

( 1.9  <  Rp  < 3.4 )  .  When  applying  the  A-coef  f  icient  to 

characteristics  of  multiple  thermohaline  staircases,  the 
average  flux  is  6.65  WmA-2  for  Kelley's  formulation  and 
3.717WmA-2  for  Marmorino  and  Caldwell's.  The  range  of 
fluxes  is  obtained  using  ITP  1  data  from  April  of  2006. 
Average  AT  ,R,  and  H  values  were  calculated  in  20  meter 

increments  from  220  to  360  meters  (except  for  320-360  meters 
where  steps  were  too  large  to  make  accurate  calculations 
from  a  20  meter  interval) .  This  range  covers  1211  steps. 
Table  3  shows  the  data  averaged  by  layer. 


Depth  in  meters 

220-240 

240-260 

260-280 

280-300 

300-320 

320-360 

AT 

0.0426 

0.0405 

0.0423 

0.0635 

0.0439 

0.0368 

AS 

0.0132 

0.0125 

0.0118 

0.0202 

0.0145 

0.0138 

Rp 

6.1626 

5.5163 

4.6419 

4.9691 

4.9435 

5.6985 

H  (meters) 

2.3520 

2.6326 

3.1538 

4.4798 

6.0994 

5.9379 

Total  Steps 

317 

306 

275 

157 

70 

86 

Kelley's  Flux  Law 

HF  w/A=1 

0.1570 

0.1637 

0.2134 

0.3327 

0.2045 

0.0126 

HFw/A=  19.6078 

3.0783 

3.2091 

4.1853 

6.5227 

4.0099 

0.2472 

HF  w/A=30 

4.7090 

4.9099 

6.2648 

9.9797 

6.1351 

4.1831 

HF  w/A=31.01522 

4.8692 

5.0761 

6.6201 

10.3175 

6.3427 

0.3910 

Marmorino  and  Caldwell’s 
Flux  Law 

HF  w/A=1 

0.1524 

0.1606 

0.2176 

0.3370 

0.2076 

0.1366 

HF  w/A=8.75 

1.3333 

1.4056 

1 .9039 

2.9488 

1.8162 

1.1953 

HF  A=15.75 

2.4000 

2.5300 

3.4270 

5.3079 

3.2691 

2.1515 

HF  w/A=1 7.955 

2.7360 

2.8842 

3.9068 

6.0510 

3.7268 

2.4527 

Table  3.  ITP  1  data  for  April  2006  in  20-40  meter  layer 

increments . 
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Table  4  uses  the  averaged  data  from  Table  1  to  provide 
a  range  of  possible  heat  fluxes  in  the  Arctic.  The  range 
varies  significantly  based  on  the  size  of  the  coefficient  A. 
This  suggests  that  further  study  is  necessary  to  determine 
the  validity  of  the  4/3  flux  law  when  applying  it  to  ocean 
conditions . 


Table  4.  Heat  Fluxes  in  the  Beaufort  Gyre  using  the 

calibrated  flux  laws. 
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V.  ACOUSTIC  APPLICATIONS 


A.  EFFECTS  OF  THERMOHALINE  STEPS  ON  ACOUSTIC  TRANSMISSION 

The  "fine-structure"  and  "microstructure" 
characteristics  of  diffusive  convection  described  above  were 
rarely  considered  to  have  an  impact  on  acoustic 
transmission.  However,  with  recent  high-resolution  ocean 
data  illustrating  the  depth  and  range  that  thermohaline 
staircases  cover,  these  small-scaled  phenomena  are  given  new 
consideration.  Chin-Bing  et  al .  (1994)  studied  the  effects 
of  the  thermohaline  staircase  step-structure  on  ocean  sound 
velocity  profiles  (SVP's)  and  then  performed  acoustic 
propagation  simulations  using  the  oceanographic  data. 

The  simulations  done  by  Chin-Bing  et  al . ,  used  data 
collected  from  the  northeast  coast  of  South  America  in  the 
fall  (October-November )  of  1985.  The  staircases  studied 
were  the  result  of  salt-fingers  and  therefore  had  different 
characteristics  (i.e.,  step  size  and  density  gradients)  from 
the  present  studies  of  thermohaline  staircases.  However, 
these  results  led  to  the  initial  consideration  of  acoustic 
propagation  through  diffusive  convective  layers. 

B.  INITIAL  APPLICATION  TO  DIFFUSIVE  CONVECTION 

Wide-angle  capable,  research-level.  Parabolic  Equation 
(PE)  models  were  used  to  simulate  long-range  transmission 
loss  through  ocean  thermohaline  steps.  The  wide-angle  models 
were  used  due  to  the  deep  ocean  waveguide  which  supports 

fully  retracted  propagation  for  angles  of  ±15°  .  Due  to  the 
sensitivity  of  PE  models  to  range  step  and  depth  step,  the 
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model  was  run  decreasing  range  and  depth  steps  until  the 
results  converged  on  a  stable  answer.  This  led  to  the  use 
of  range  and  depth  step  sizes  of  1/24  the  wavelength.  The 
frequency  choice  of  50Hz  gave  an  acoustic  wavelength  on  the 
same  order  as  the  step-size  in  the  SVP .  The  source  depth 
was  400  meters,  placing  it  in  the  mid-depth  of  the  steps. 
Results  were  plotted  for  receiver  depths  of  200,  400,  and 

800  meters,  that  is,  above,  in,  and  below  the  step  region 
respectively.  In  order  to  determine  the  importance  of 
transmission  loss  through  the  staircases,  a  comparison  was 
made  between  the  transmission  loss  of  the  high-resolution 
step  sound  velocity  profiles  (SVP)  and  lower  resolution  data 
only  taken  at  standard  depths . 


WIEN  UAL  DENSITY. 
26V?_2GJL27.1 

salinity,  s 

342  150  3S2  3S.4  15.6  SS.B  36.0 
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POTENTIAL  TEMPERATURE,  # 


b) 


Comprison  of  thrw  SVP’s  from  tlrermohaline  step 
standard  depth  data,  and  historical  data 


Figure  29.  a)  Typical  temperature  and  salinity  plot  for 

salt  finger  region.  Salinity,  Potential  Temperature 
and  Density  plots  show  thermohaline  step 
characteristics.  b)  Sound  velocity  profiles  showing 
the  step  data,  standard  depth  profiles  and  historical 
averaged  sound  velocity  profile  data.  From  Chin-Bing  et 

al.  (1994)  . 
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The  resulting  simulation  showed  that  when  the  receiver 
was  above  the  step  region  there  was  very  little  difference 
when  comparing  the  transmission  loss  of  the  two  SVP  data 
sets.  The  same  determination  was  made  when  the  receiver  was 
under  the  staircase  structure.  There  was,  however,  a 
significant  difference  within  the  step  region--5  decibels 
over  a  10-kilometer  range.  The  greatest  impact  on  the  ray 
paths  were  when  both  the  source  and  receiver  were  at  400 
meters,  both  were  within  the  staircase  structure. 

The  numerical  experiment  then  went  on  to  test  the 
effects  of  the  staircases  on  backscatter.  Using  the  same 
50Hz  source,  a  source  depth  of  100  meters  below  the  surface 
was  chosen.  Chin-Bing  et  al .  found  that  in  the  low 
frequency  regime  of  the  study,  the  staircases  tend  to  have  a 
negligible  effect;  however,  they  felt  as  the  frequency  of 
the  acoustic  field  increased  so  should  the  magnitude  of  the 
backscatter.  They  surmise  that  given  frequencies  in  the 
sonar  range  i.e.  3.0  kHz,  the  backscatter  from  thermohaline 
steps  could  be  significantly  important. 

C.  THE  ARCTIC  CASE 

From  the  above  research,  we  thought  it  prudent  to  test 
the  impact  of  the  diffusive  convection  form  of  thermohaline 
staircases,  given  the  probability  of  increased  interest  in 
the  Arctic.  The  raw  data  from  ITP  profile  1013,  chosen  for 
its  well  formed  staircase  structure,  was  used  for  basic 
experimentation  (Figure  30)  . 
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Figure  30.  Temperature  and  depth  profile  used  for 

acoustic  simulation . 

Raw  temperature  data  and  converted  salinity  (from 
conductivity)  data  was  used  for  the  idealized  experiment. 
Using  a  basic  fit  of  the  data,  the  staircase  structure  was 
smoothed  from  the  profile  in  order  to  offer  a  comparison  in 
the  effects  of  the  microstructure  on  sound  propagation. 
Using  Del  Grosso's  Equation  (Del  Grosso  1974),  sound  speed 
profiles  were  calculated  for  simulation  purposes.  Figure  31 
shows  the  two  sound  velocity  profiles.  The  only  region  of 
the  profiles  that  was  smoothed  was  the  staircase  region,  so 
that  effects  of  other  factors  were  not  included  in  the 
calculated  impact . 
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S  S  PI  01 3  /  smoothed  /with  and  without  staircases 


Figure  31.  Sound  Velocity  Profile  from  ITP  1  rawl013 

data.  Red  profile  smoothes  staircase  structure,  blue 
includes  staircases.  The  graphic  on  the  right  is  an 
enhancement  to  the  SVP  so  that  the  staircase  region  is 

more  visible. 

The  parameters  for  the  simulation  were  a  planar, 
perfectly  reflecting  surface.  (This  is  not  the  case  for 
under  the  ice  sheet,  where  there  would  be  increased 
scattering  of  sound  as  the  ray  paths  encountered  the  ice.) 
The  bottom  was  simulated  as  a  17  meter  thick  sediment  layer 
with  a  bottom  half  space.  A  sandy  bottom  type  was  assumed. 
Due  to  the  upward  refracting  characteristics  if  the  Arctic, 
the  bottom  type  did  not  play  a  role  in  propagation. 

The  BELLHOP  MODEL,  which  is  a  Gaussian  beam  propagation 
model  using  coherent  addition  of  ray  paths,  was  used  in 
order  to  keep  the  model  run  time  reasonable.  The  3.5  kHz 
source  was  placed  at  a  depth  of  ten  meters.  Figures  32  and 
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Depth  (m) 


33  are  the  resulting  ray  traces  with  transmission  loss  (TL) 
referenced  dB  re  1  micro  Pascal. 
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Figure  32 . 


Sound  propagation  through  SVP  with  no 
staircases . 
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Figure  33.  Sound  Propagation  model  through  SVP  with 

staircase  structure  included. 


Logarithmic  scaling  makes  it  difficult  to  easily 
compare  the  two  simulations,  so  the  difference  between  the 
two  cases  is  plotted  to  provide  clarity.  From  Figure  34,  it 
can  be  determined  that  within  a  range  of  10000  meters, 
(range  of  operational  detection) ,  there  is  a  significant 
impact  from  the  staircases,  with  differences  in  TL  up  to  15 
dB .  It  can  be  noted  that  there  is  very  little  difference  in 
the  ray  traces  within  the  first  4000  meters.  This  is  an 
expected  result  because  in  order  for  the  staircases  to  have 
an  impact,  the  sound  has  to  be  refracted  so  that  it 
encounters  the  staircases  at  near-horizontal  angles  (a 
reasonable  assumption  for  any  source  placement  above  the 
staircases  because  of  upward  refraction)  The  difference  in 
TL  using  this  idealized  model  is  between  -15  and  15  dB. 
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Figure  34.  Difference  in  TL  between  profiles  with  and 

without  staircases . 

Considering  the  idealized  situation,  it  is  reasonable 
to  assume  that  scattering  from  the  ice,  or  other  ocean 
sounds  speed  structure,  could  reduce  or  enhance  the  range  of 
impact  of  the  staircases.  It  should  be  noted  that  even  a  3dB 
difference  can  cause  a  significant  impact  on  detection  and 
should  be  considered.  Further  research  is  needed  to  realize 
that  full  impact  of  these  staircases  on  sound  propagation. 
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VI .  CONCLUSIONS 


Among  various  aspects  of  the  diffusive  convection 
explored  in  this  study,  the  most  significant  was  the 
development  of  the  layer-merging  theory  for  thermohaline 
staircases.  The  theory  based  analysis  of  the  Beaufort  Gyre 
made  it  possible  to  evaluate  existing  flux  laws  for  the 
ocean.  From  ice-tethered  profiler  data,  there  is  evidence 
that  merging  follows  Radko  (2007)  B-merger  theory.  Applying 
the  theory  directly  to  in  situ  data,  it  became  evident  that 
the  pattern  of  the  normalized  growth-rate  was  the  same. 
However,  it  was  determined  that  the  laboratory  derived  flux 
laws  cannot  be  directly  applied  to  Arctic  staircases.  A 
coefficient  is  necessary  to  convert  from  laboratory  to  the 
ocean.  The  coefficient  is  determined  by  analyzing  the 
merging  pattern  of  diffusive  layers,  and  relating  it  to 
theory  based  flux  laws.  The  proposed  theoretical  model  of 
layer  merging  events  is  consistent  with  the  field  data. 

Based  on  this  consistency,  a  coefficient  "A"  0(1)  can 
be  determined  and  applied  to  better  correlate  laboratory 
flux  laws  with  ocean  data.  Using  different  expressions  for 
the  R  dependent  coefficient  of  the  4/3  flux  law  (Kelley, 

1990;  Marmorino,  1976),  we  estimate  that  the  average  flux 
from  the  thermohaline  staircases  is  from  1-6  WmA-2.  This 
range  of  fluxes  is  significantly  higher  than  if  the  original 
4/3  flux  law  is  applied.  Inconsistencies  in  the  results  of 
laboratory-derived  fluxes  provide  significant  room  for  error 
in  direct  extrapolation  to  the  ocean.  Therefore,  the  higher 
fluxes  derived  through  merging  theory  cannot  be  discounted. 
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Fluxes  in  the  Arctic  are  currently  of  great  concern.  As 
Turner  and  Veronis  (2004)  suggest  through  their  laboratory 
experimentation  on  ice  melt,  double-diffusive  processes  may 
be  somewhat  responsible  for  melting  Arctic  ice.  While  it 
has  been  suggested  that  fluxes  in  the  Beaufort  Gyre 
staircases  cannot  penetrate  the  halocline,  it  is  entirely 
possible  that  regions  of  double-diffusive  convection,  where 
the  halocline  is  weaker,  could  have  a  greater  effect.  It  is 
also  noted  that  the  heat  budget  for  the  Arctic  has 
approximately  1/3  of  its  fluxes  not  accounted  for  by  solar 
input  (Perovich  et  al .  2001)  .  While  evidence  is 
inconclusive  as  to  the  exact  contribution  of  convective 
diffusive  fluxes  to  mixed  layer  over-all  flux,  it  would  be 
negligent  not  to  consider  the  possibility  of  its 
contribution . 

In  the  cycle  of  global  climate  change,  attention  is 
once  again  being  diverted  to  the  Arctic.  The  accessibility 
of  resources  and  shipping  routes  that  a  melted  ice  cap  would 
provide  necessitate  an  operational  focus  on  the  impact  of 
microstructures  on  sound  propagation.  An  idealized 
propagation  model  was  used  to  determine  whether  or  not 
further  research  was  necessary,  and  found  that  under  normal 
operational  conditions,  there  was  an  impact  from  the 
staircase  structure  of  up  to  30dB.  While  further  research 
could  increase  or  decrease  this  magnitude,  the  difference  in 
transmission  loss  should  be  considered  in  future  Arctic 
operations,  and  should  be  further  studied  to  determine  the 
actual  magnitude  of  impact  from  this  structure. 
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